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Stem cell fate and function are dynamically modulated by the interdependent relationships between 
biochemical and biophysical signals constituting the local 3D microenvironment. While approaches to 
recapitulate the stem cell niche have been explored for directing stem cell differentiation, a quantitative 
relationship between embryonic stem cell (ESC) morphogenesis and intrinsic biophysical cues within 
three-dimensional microtissues has not been established. In this study, we demonstrate that mesenchymal 
embryonic microtissues induced by BMP4 exhibited increased stiffness and viscosity accompanying 
differentiation, with cytoskeletal tension significantly contributing to multicellular stiffness. Perturbation 
of the cytoskeleton during ESC differentiation led to modulation of the biomechanical and gene expression 
profiles, with the resulting cell phenotype and biophysical properties being highly correlated by multivariate 
analyses. Together, this study elucidates the dynamics of biophysical and biochemical signatures within 
embryonic microenvironments, with broad implications for monitoring tissue dynamics, modeling 
pathophysiological and embryonic morphogenesis and directing stem cell patterning and differentiation. 



The biochemical composition and physical structure of tissue microenvironments are critical regulators of 
stem cell fate and function 1,2 . In vivo, native stem cell niches maintain the self-renewing state through short 
range soluble signaling 3 and extracellular signals 4 , such as direct anchorage via intercellular and/or matrix 
adhesions. Conversely, dysregulation of niche elements leads to differentiation and migration of resident stem cell 
populations 5 . Therefore, the synergy of biophysical and biochemical signals and the complexity of the local 
microenvironment have been increasingly recognized as critical mediators of the delicate balance between 
self-renewal and differentiation. However, while engineering approaches often aim to dissect the relative influ- 
ence of individual stimuli, the signals regulating stem cell differentiation are inherently interrelated, motivating 
the multiparametric analysis of stem cell structure and phenotype accompanying dynamic cell fate decisions. 

The pre-implantation embryo is one of the most dynamic tissue microenvironments, whereby the patterning of 
structurally and functionally distinct tissue structures, such as liver, skin, and heart arise from physically adjacent, 
yet precise spatially defined, regions during development 6 . In particular, embryonic gastrulation mediates the first 
somatic cell fate decisions, and results in spatially distinct localization of the three germ lineages - endoderm, 
ectoderm, and mesoderm. During gastrulation, the cells of the epiblast that will comprise the mesoderm lineage 
undergo a rapid and concerted morphogenesis 7,8 through which the tightly adherent, epithelial cells become less 
adhesive mesenchymal cells, in a process known as the epithelial-to-mesenchymal transition (EMT) 9 . The 
dynamics and patterning of EMT are tightly regulated by the biophysical and biochemical tissue microenviron- 
ment, through signals that maintain or abrogate the apical-basal polarity of epithelial cells. While the role of 
biochemical signals are well established in the induction of EMT, the resulting changes in cell shape that occur 
during gastrulation are a fundamentally biomechanical process, through which forces arising as a result of 
adhesion remodeling are transmitted via changes in cytoskeletal tension 1012 . Biophysical signals also act in 
conjunction with biochemical cues to mediate different stages of embryonic morphogenesis, including mesoderm 
intercalation 13 , dorsal elongation 14 , and neural tube closure 15 . Moreover, mechanical cues alter the responsiveness 
of cells to chemical signals in order to mediate proliferation and dictate organ size 16 , thereby demonstrating the 
complexity of deciphering signaling mechanisms regulating tissue patterning. 
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While the phenotypic 17 and mechanical 18 changes arising during 
embryonic development have been characterized through ortho- 
gonal approaches in various species 6,19,20 , several questions remain 
regarding the role of cellular biomechanics in determining cell fate, 
particularly during mammalian embryogenesis. The differentiation 
of pluripotent embryonic stem cells (ESCs), the in vitro analog to the 
inner cell mass of pre-implantation embryos, parallels many of the in 
vivo cell fate transitions, reflecting the utility of ESCs as a model for 
systematically studying mammalian embryonic development 21 . ESC 
differentiation approaches have aimed to recapitulate the biochem- 
ical and/or biophysical milieu of the gastrulating embryo, including 
the perturbation of developmentally relevant signaling pathways 22 
and the manipulation of the composition and mechanical properties 
of adherent substrates 23,24 . However, when assembled as three- 
dimensional aggregates, ESCs undergo morphogenic events, includ- 
ing the deposition of extracellular mediators of EMT, such as 
hyaluronan and versican 25 , which initiates dynamic, development- 
ally relevant processes, thereby offering unique opportunities to 
study biomechanics in parallel with changes in pluripotent cell fate 
and morphogenesis. 

The objective of this study, therefore, was to define the intrinsic 
cellular biophysical characteristics and phenotypic changes that arise 
simultaneously during cell fate specification and morphogenesis of 
three-dimensional pluripotent stem cell microenvironments. A mul- 
tifaceted approach was employed to assess the dynamics of morpho- 
genesis and differentiation, characterize the biomechanical 
properties of 3D microtissues, and determine correlative associations 
between gene expression and mechanical signatures via multivariate 
data modeling. Overall, this study demonstrates that morphogenesis 
arises via EMT during mesoderm differentiation of pluripotent 
microtissues, resulting in dynamic temporal changes in phenotype 



that are highly correlated with biomechanical properties. Ultimately, 
understanding the biophysical changes accompanying epithelial-to- 
mesenchymal transition of ESCs has broad implications to reveal 
physical mechanisms underlying mammalian embryonic develop- 
ment, tissue homoeostasis, pathological remodeling and cancer 
metastasis. 

Results 

BMP4 induces mesoderm differentiation in ESC microtissues. 

Homogeneous populations of EBs were formed via forced aggre- 
gation and maintained in controlled hydrodynamic suspension 
culture, as described previously 26 , in order to systematically examine 
three-dimensional stem cell morphogenesis. As expected, supple- 
mentation with BMP4 (10 ng/mL) altered the EB morphology 
(Fig. la), as evidenced by changes in the EB circularity and optical 
density (arrows), and overall gene expression profile (Fig. lb), 
resulting in the increased expression of genes related to mesoderm 
lineages (Supplementary Table SI) at days 4 and 7 of differentiation 
compared to those maintained in basal, serum-free culture conditions. 
Although decreases in the pluripotency markers Sox2 and Oct-4 did 
not differ significantiy between the two culture conditions, markers of 
hematopoietic mesoderm {Gata.2, Hba-x) were increased by BMP4 
treatment (Gata2 p = 0.07; Hba-x p = 0.03), and genes related to 
endoderm (Foxa2) and ectoderm (Fgf5) were significantly decreased 
(Foxa2 p = 0.04; Fgf5 p = 0.004), thereby illustrating the mor- 
phogenic influence of BMP4 in directing mesoderm differentiation 
in 3D cultures, analogous to previous reports using monolayer diffe- 
rentiation techniques 27,28 . 

ESCs undergo mesenchymal morphogenesis during 3D 
mesoderm differentiation. Consistent with the gross changes in 
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Figure 1 | Induction of mesoderm differentiation within EBs via BMP4 treatment, (a) Phase contrast images demonstrate similar gross morphology of 
EBs cultured in basal, serum-free media or media supplemented with BMP4 after 4 days of differentiation, with differences in circularity (basal) and 
optical density (BMP4), particularly apparent after 14 days of differentiation (arrows), (b) Gene expression profiles also differed between EBs cultured in 
different soluble culture environments after 4 and 7 days of differentiation, (c) resulting in similar levels of pluripotency factors, with divergent expression 
of markers for hematopoietic mesoderm (Gata2, Hba-x), endoderm (Foxa2), and ectoderm {Fgf5). Scale bar = 200 urn; n = 3; * = p < 0.05 using two- 
way ANOVA and post-hoc Tukey analysis. 
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Figure 2 | Mesenchymal morphogenesis of ESCs within EBs. (a) EBs exhibited dynamic remodeling of the microenvironment through 14 days of 
differentiation, with distinct differences apparent between EBs (H&E staining) in different soluble media formulations, (b) After 7 days of differentiation, 
mesenchymal-like regions (arrows) within EBs exhibited positive staining for alcian blue (blue) and a-SMA (brown/black), whereas safranin-O (red/ 
brown) stained positively within epithelial-like regions and was largely absent within the mesenchymal-like regions of BMP4-treated EBs. Scale bars (A) 
= 200 p,m, (B) = 50 (im; n = 3; Nuclei counterstained with: fast green (safranin-O), nuclear fast red (alcian blue), and hematoxylin (a-SMA). 



EB morphology, distinct differences in cellular morphology and 
organization within EBs were observed during the course of 
differentiation, in a soluble factor- dependent manner (Fig. 2a). 
While EBs within both conditions exhibited similarly homogenous 
populations of cells with epithelial morphologies initially (day 4), 
biochemical-mediated differentiation led to changes in cellular 
organization and morphogenesis, with basal EBs continuing to 
exhibit tightly packed, epithelial cell clusters, in contrast to BMP4- 
treated EBs, in which populations of elongated, mesenchymal-like 
cells were apparent at day 7 of differentiation. Interestingly, the 
mesenchymal morphology observed at day 7 was no longer 
apparent by 14 days, indicating a transition phase accompanying 
the specific stage of mesoderm differentiation. Consistent with the 
GAG-rich pericellular matrix of undifferentiated ESCs, the EBs 
maintained in basal conditions exhibited positive safranin-O 
staining. In addition to the morphological changes, BMP4-treated 
EBs exhibited a glycosaminoglycan (GAG) rich matrix at day 7, with 
spatial variation in composition, as evidenced by positive alcian blue 
(pH 2.5) and negative safranin-O staining within the interstitial 
space between the regions of mesenchymal cells (Fig. 2b). The 
putative mesenchymal cells also expressed a-smoofh muscle actin 
(a-SMA), an established marker of EMT 29 . Taken together, EB 
mesoderm differentiation induced by BMP4 treatment results in 
EMT morphogenesis to yield characteristic mesenchymal cell 
populations. 

Biophysical microtissue properties are dynamically modulated 
during ESC differentiation. To characterize the physical 
remodeling of the EB microenvironment during ESC differentiation, 
biomechanical microtissue properties were measured via micron-scale 
parallel plate compression (Supplementary Fig. Sla). Descriptive 
physical characteristics, including moduli, were calculated by 
measuring the viscoelastic creep (change in deformation under 
constant force), which was fitted to a linear viscoelastic model 
(average R 2 = 0.97; Supplementary Fig. Slb,c). EBs treated with 
BMP4 exhibited increased viscous resistance to deformation (p = 
0.0003), as evidenced by decreased creep compared to those 



differentiated in basal conditions at day 7 of differentiation (Fig. 3a). 
BMP4-treated EBs were 1.5-fold stiffer than basal EBs after 7 days of 
differentiation (p = 0.004; Fig. 3b). Therefore, morphologically 
distinct EBs also exhibited biomechanical differences, both in terms 
of viscoelasticity and microtissue stiffness. Interestingly, the size 
distribution of EBs was altered upon BMP4 treatment, which led to 
increased dispersity after 7 days of differentiation (median absolute 
deviation, MAD = 31 urn for basal EBs; MAD = 62 um for BMP4 
EBs). While the stiffness and size of basal EBs were not strongly 
correlated (Fig. 3c), those differentiated with BMP4 demonstrated 
an inverse relationship (p < 0.001), with smaller BMP4-treated EBs 
exhibiting increased stiffness (Fig. 3d), thereby suggesting that smaller 
EBs may undergo more robust mesenchymal responses to BMP4, 
potentially due to transport-related phenomena 30 . In addition to the 
increased stiffness of the BMP4-treated population of EBs, the 
correlation between aggregate size and stiffness reflects the physical 
relationship between morphogenesis and biomechanical properties, 
thereby providing a novel quantitative metric to assess the state of 
EB differentiation. 

Moreover, the biomechanical properties of EBs were dynamically 
modulated during the course of differentiation, whereby significant 
differences in the stiffness and viscosity were measured over time, 
both within and between individual biochemical differentiation con- 
ditions (Fig. 3e,f and Supplementary Fig. S2). Specifically, basal con- 
ditions exhibited decreased stiffness between days 2 and 4 of 
differentiation (p = 0.04), whereas BMP4-treated EBs exhibited sim- 
ilar stiffness throughout the first 7 days of differentiation; the stiff- 
ness of basal and BMP4-treated EBs was significantly increased at 
day 14 compared to day 4 (basal p = 0.01; BMP4 p < 0.001), thus 
indicating an overall increase in stiffness during differentiation. Both 
conditions also exhibited similar initial viscosity dynamics, whereby 
the viscous response decreased (i.e. increasing deformation time, or 
creep time constant) during the initial 7 days of differentiation (basal 
p = 0.04; BMP4 p = 0.04); however, the viscosity of EBs from basal 
conditions continued to decrease through day 14 of differentiation, 
whereas the viscosity of BMP4-treated EBs increased between days 7 
and 14 of differentiation, resulting in significantly increased viscosity 



SCIENTIFIC REPORTS | 4:4290 | DOI: 1 0. 1 038/srep04290 



3 



100 
80 
60 



S. 40 



20 
0 

120 



a. 



■5 80 

O 

E 



40 



CL 

cl 
< 



Creep deformation 



t 




p = 0.0003 



T 



Basal 



+BMP4 



Relaxed modulus, 



9p = 0.004 



T 



Basal 



+BMP4 



C 120 



CO 
CL 



T3 
O 

E 



80 



* 40 



co 
a. 

Cl 
< 



0 



Basal 



to o 



o 



o o(j o 

o c? r = -0.25 
(N.S.) 



— 40 
in 



30 



50 ., Creep time constant, T a 

-o - Basal 

- BMP4 T 

A 



200 300 400 500 600 



20 



120 



i/l 



■5 80 

O 

E 



Cl 

CL 
< 



40 



0 



BMP4 
• • • 

• » .* 

r = -0.79 • 
(p = 6x10" 7 ) 



f 






300 


"co" 




</> 




_= 
3 


200 


T3 




O 




£ 




<u 


100 


1_ 

to 

CL 




CL 




< 






0 



0 2 4 6 8 10 12 14 



Relaxed modulus, E M 



200 300 400 500 600 

EB diameter (pm) 



f -£---o-' 



,4 



.---5 



0 2 4 6 8 10 12 14 

Day of Differentiation 



Figure 3 | Dynamic modulation of EB mechanical properties. EBs exhibited differences in (a) viscoelastic creep and (b) modulus after 7 days of 
differentiation, including (c,d) a significant correlation between modulus and EB diameter within BMP4-treated EBs. In addition, (e) modulus and 
(f) viscoelastic time constant characteristics were modulated as a function of time, with differences between culture conditions after 14 days of 
differentiation. (a,b) gray shaded regions represent the upper and lower quartiles while the red center line denotes the median of the populations, 
(a-d) n = 27 EBs from 5 independent experiments; (e), (f) n = 6 EBs; * = p < 0.05 using two-way ANOVA and post-hoc Tukey analysis. 



of BMP4-treated EBs compared to basal EBs (p < 0.001). Overall, 
the biomechanical properties of multicellular ESC aggregates change 
dynamically during differentiation, with BMP4-treated EBs stiffen- 
ing and exhibiting increased viscosity as cells commit toward meso- 
dermal lineages and undergo EMT. 

Cytoskeletal tension regulates biophysical microtissue environments. 

The actomyosin cytoskeleton contributes to the primary structural 
mechanics of cells and is particularly important in mediating cell 
shape, force generation, and stiffness in embryonic tissues 31 . There- 
fore, we hypothesized that the biophysical and phenotypic changes 
arising during three-dimensional pluripotent stem cell morpho- 
genesis and differentiation were mediated by cytoskeletal tension. 
The actin filament organization of cells comprising basal serum- 
free EBs retained a largely cortical structure at all stages of diffe- 
rentiation, consistent with an epithelial cell phenotype; however, 
BMP4-treated EBs exhibited decreased assembly of F-actin, with 
punctuated cortical structures after 4 days of differentiation (Supple- 
mentary Fig. S3). Moreover, consistent with the mesenchymal 
morphology, cells within BMP4-treated EBs exhibited characteristic 
elongated cellular morphologies and filamentous cytoskeletal archi- 
tecture by day 7 of differentiation, which implicates the cytoskeleton 
in mediating microtissue stiffness during EB morphogenesis. Disrup- 
tion of the Rho associated protein kinase (ROCK) pathway with a 
small molecule inhibitor (Y27632; 10 uM; 1 hour) 32 (Fig. 4a) acutely 
decreased the stiffness of EBs from both differentiation conditions at 
all time points examined (p < 0.03 for all conditions), indicating that 
the actomyosin cytoskeleton significantly contributes to the stiffness 
of EBs (—50%) (Fig. 4 b,c). In addition, despite similar EB stiffness, 
the ROCK-mediated contractility of BMP4- treated was increased 
compared to basal EBs after 2 days of differentiation (p < 0.001), 
demonstrating the emergence of unique biophysical traits at early 
stages of differentiation, even before morphological or phenotypic 
changes could be observed (Fig. 4d). Moreover, although the 



magnitude of stiffness increased during differentiation, the relative 
contribution of cytoskeletal tension to the overall modulus of BMP4- 
treated EBs decreased, indicating an increasing role for other 
structural mechanisms (e.g. nuclear, extracellular) at later stages of 
mesenchymal morphogenesis. Altogether, cytoskeletal tension signi- 
ficantly contributes to ESC microtissue mechanics and the actomyo- 
sin cytoskeleton is remodeled during mesenchymal morphogenesis. 

While the cytoskeletal structure directly mediates tissue stiffness, 
the dynamics of actin polymerization have also been implicated in 
morphogenic processes, such as EMT 33 . Therefore, the actin polymer- 
ization was perturbed via treatment with an agonist (Jasplakinoide, 
Jas; 50 nM) 34 or antagonist (latrunculin B, LatB; 100 nM) 35 in addi- 
tion to perturbation of ROCK-mediated contractility with Y27632 
during the initial period of ESC morphogenesis (days 4-7) in order 
to assess the impact on biomechanics, morphogenesis, and cell 
phenotype (Fig. 4e). BMP4- treated EBs exhibited decreased viscosity 
in response to Jas and Y27632, compared to untreated (Jas p < 0.001; 
Y27632 p = 0.01) EBs (Fig. 4f); in contrast, the viscous response of 
basal EBs remained unchanged across all conditions. Despite opposite 
mechanisms of action, the decreased viscosity in response to both Jas 
(agonist) and Y27632 (antagonist) demonstrates additional time- 
dependent biophysical properties, which may arise from feedback 
or compensatory changes arising from dynamic changes in cellular 
organization and phenotype, as well as remodeling of adhesions, or 
modulation of the spatial distribution or relative contributions of 
intracellular structural elements, including actin, myosin, and inter- 
mediate filaments, in response to cytoskeletal perturbations. The kin- 
etics of the viscous response, reflected by the creep time constant, 
were increased by treatment with Y27632 in BMP -4 and basal con- 
ditions (Fig. 4g), suggesting that inhibiting ROCK-mediated contract- 
ility decreased aggregate viscosity. In addition, Jas increased the 
stiffness of EBs cultured in basal, serum-free media compared to 
untreated EBs, consistent with its known role stimulating actin poly- 
merization (p < 0.001); however, the stiffness of BMP4-treated EBs 
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Figure 4 | Influence of cytoskeletal tension on EB mechanical properties, (a) Treatment of EBs with the ROCK inhibitor Y27632 for 1 hour significantly 
decreased the modulus of (b) basal and (c) BMP4-treated EBs and (d) demonstrated differences in the relative contribution to the overall naive EB 
modulus at 2 and 7 days of differentiation, (e) In addition, EBs were differentiated in the presence of cytoskeletal agonists and antagonists, and exhibited 
changes in the (f) viscoelastic creep and (g) time constant responses, as well as (h) modulus after 7 days of differentiation, n = 6 EBs; * = p < 0.05 
compared to basal; $ = p < 0.05 compared to untreated, using two-way ANOVA and post-hoc Tukey analysis. 



was not affected by Jas (Fig. 4h). In contrast, Y27632 significantly 
decreased the stiffness of BMP4 EBs compared to untreated condi- 
tions (p < 0.001), but did not affect EBs differentiated in basal con- 
ditions. Interestingly, although LatB has been reported to decrease 
actin polymerization, BMP4-treated EBs were significantly stiffer 
when cultured in LatB compared to untreated EBs (p < 0.001). 
Overall, these studies indicate that the biomechanical properties of 
pluripotent microtissues undergoing differentiation and morphogen- 
esis are significantly altered by the dynamics of actin polymerization 
contributing to cytoskeletal tension. 

Biomechanics predict microtissue phenotypic profiles. The impact 
of biomechanical forces on stem cell phenotype has been increasingly 
appreciated as an exogenous means to direct cell fate 24,36,37 ; however, 
a quantitative relationship between ESC phenotype and biophy- 
sical properties of three-dimensional microtissues has not been 
established. Therefore, the morphology and phenotype of ESCs 
were assessed in parallel with the changes in microtissue stiffness 
upon treatment with small molecules. The cytoskeletal agonist and 
antagonists did not induce marked morphological changes, based on 
gross histology (Supplementary Fig. S4); however, subtle differences 
were noted, as LatB-treated BMP4 EBs exhibited more epithelial-like 
cells, similar to basal EBs (Supplementary Fig. S4s), whereas Jas and 
Y27632 treatment induced more mesenchymal morphogenesis 
within BMP4 EBs (Supplementary Fig. S4q,r). In addition, treat- 
ment with the small molecules altered the gene expression profiles 
of both basal and BMP4 EB populations after 7 days of differen- 
tiation; hierarchical clustering indicated marked phenotypic 
differences for Jas treated EBs from both conditions relative to 
untreated EBs (Fig. 5a). Similarly, a PLS model (R 2 Y = 0.70, Q 2 = 



0.10; 3 significant PCs) illustrated the variance in gene expression 
signatures due to BMP4 treatment and cytoskeletal perturbations, 
which were distinctly separated along the first and second principal 
components, respectively (Fig. 5b); although the biochemical 
morphogen treatment in the presence or absence of BMP4 
dominated the phenotypic gene expression, subtle variations in the 
phenotypic profile were induced by cytoskeletal perturbations, 
thereby indicating the utility of PLS modeling for capturing 
variations not detectable through traditional analytical methods, 
such as hierarchical clustering. Consistent with the histological 
observations, LatB-treated EBs exhibited phenotypic profiles 
indicative of decreased EMT, similar to basal EBs, whereas Jas- 
treated BMP4 EBs exhibited increased mesoderm differentiation 
and mesenchymal morphogenesis, as evidenced by the position 
along the first principal component. Increased expression of 
pluripotency genes such as Nanog, Sox-2, and Oct-4, indicated that 
LatB treatment delayed the differentiation kinetics of ESCs for both 
culture media conditions. Interestingly, the same mesendoderm 
genes (Hbb-y, Soxl7, Nkx2.5, Fgf5) increased by Jas treatment of 
BMP4 EBs were decreased by Jas in basal EBs, thus demonstrating 
that actin polymerization was not sufficient to induce a mesendo- 
derm gene expression profile within basal EBs (Fig. 5b). 

Based on the biochemical-mediated variability between experi- 
mental groups captured by the first principal component, separate 
PLS models were constructed based upon the gene expression (input, 
X) and mechanical properties (output, Y) to distinguish the cellular 
responses to cytoskeletal perturbations for EBs differentiated in basal, 
serum-free media (R 2 Y = 0.70, Q 2 = 0.10; 1 significant PC) and in 
the presence of BMP4 (R 2 Y = 0.65, Q 2 = 0.30; 1 significant PC) 
(Fig. 5c). The models reflect changes in biomechanical properties as 
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Figure 5 | Prediction of microtissue stiffness and ESC phenotypes by multivariate modeling, (a) PCR array analysis of gene expression suggested 
differences in cell fate after 7 days of differentiation due to small molecule-mediated perturbation of the cytoskeleton from days 4-7. (b) In addition, the 
score plot demonstrates distinctions between culture conditions across the first and second principal components, based upon gene expression profiles, 
shown in the loadings plot, (c) Separate PLS with gene expression as inputs and mechanical properties as outputs demonstrated a high degree of 
correlation between mechanical properties and phenotypic traits, (d) In addition, PLS models based upon mechanical characteristics demonstrated 
different degrees of correlation for the associated gene expression profiles when trained with individual or full biomechanical profiles. 
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a function of the gene expression profile, thus enabling a quantitative 
relationship between biophysical microtissue environments and ESC 
differentiated phenotypes. Overall, the gene expression profile was 
highly correlated with several of the mechanical characteristics; how- 
ever, the strength of correlations varied between individual biome- 
chanical parameters and across differentiation conditions. The 
model theoretically predicted the EB stiffness, with R 2 coefficients 
of 0.90 and 0.75 in the basal and BMP4 conditions, respectively. In 
contrast, the viscous creep response of EBs did not correlate as highly 
with the phenotype of EBs in either condition, indicating that the 
changes in individual mechanical properties, such as viscosity, could 
not be fully explained through phenotypic changes and may also be 
influenced by extracellular factors, such as ECM deposition. To illus- 
trate the complex relationship between biophysical characteristics 
and ESC phenotype, the model was instead constructed using 
the biomechanical parameters to predict gene expression profiles 
(Fig. 5d). Overall, the biomechanical profile significantiy (back- 
ground Q 2 Y < 0.05) predicted the responses of 61.5 and 65.4% of 
the genes, respectively, which established subsets of genes responsive 
to biophysical perturbations. Moreover, the upper quartile of genes 
(Basal: Pax-6, Oct-4, Bmp4, Hba-x; BMP4: Hba-x, Hbb-y, Bmp4, 
Gata4) within each condition were predicted by the models with 
> 70% accuracy. Interestingly, Bmp4 and Hba-x, which were highly 
predicted by both models, exhibited divergent responses to small 
molecule treatments within basal- and BMP4- treated EBs, thereby 
demonstrating the capacity for multivariate analytics to identify 
putative markers highly correlated with biomechanics across differ- 
ent environmental conditions. In addition, the genes most highly 
predicted in basal EBs included those that exhibited decreased 
(Hba-x, Hbb-y), increased (Gata4), and unchanged (Bmp4) express- 
ion patterns in response to cytoskeletal perturbations, indicating that 
the mechanical model captured an array of different phenotypic 
responses. In contrast, while the mechanical model largely recapitu- 
lated the phenotypic profile of BMP4- treated EBs, several individual 
genes were not as highly correlated, which may be due to diverse 
responses arising within less uniform microtissues. The data together 
indicate that distinct phenotype profiles could be accurately 
described by multivariate analyses based upon multiple viscoelastic 
parameters. In summary, these results collectively suggest that ESC 
differentiated phenotypes are directly related through complex, mul- 
tiparametric changes in the stiffness and viscoelastic characteristics 
of EBs. 

Discussion 

This study elucidates intrinsic biophysical mechanisms mediating 
three-dimensional pluripotent morphogenesis, which can enable 
the development of novel approaches to direct the differentiation 
and patterning of functional ESC-derived microtissues for modeling 
mammalian pathophysiology and embryogenesis 38 , pharmacological 
screening 39 , and facilitating in vivo tissue repair 40 . We demonstrate 
that BMP4-mediated mesoderm differentiation of ESCs recapitulates 
aspects of embryonic development, particularly in terms of morpho- 
genesis leading to regions of mesenchymal-like cells, analogous to 
primitive streak formation during gastrulation 7,41 . Moreover, pheno- 
typic and morphogenic changes were accompanied by remodeling of 
cytoskeletal elements and modulation of biophysical properties dur- 
ing three-dimensional ESC differentiation, with the phenotypic and 
mechanical traits correlated by multivariate modeling. 

The improved understanding of three-dimensional ESC biophys- 
ical and phenotypic dynamics establishes distinctions between 
embryonic and adult tissue models and provides important princi- 
ples for the design of relevant microenvironments to direct early cell 
fate decisions. The biophysical observations noted throughout this 
study are consistent with the dramatic differences in tissue structure 
and cellular composition between adult- and embryonic-derived tis- 
sues 17 . Several previous studies have engineered the mechanical 



properties of stem cell environments through the development of 
substrates with elasticity comparable to the range of adult tissues 
(1 kPa-100 kPa) 24 ' 42 , which is 10 4 Pa less than rigid standard poly- 
styrene culture dishes. In contrast, the stiffness of morphogenic 3D 
ESC multicellular aggregates is substantially less than that of native 
tissues (Eoo < 0.1 kPa). Similar to embryonic tissues 43 , ESCs undergo 
rapid growth coincident with morphogenesis, leading to high nuclear 
density of EBs evidenced by histological analysis (Fig. 2a,b), which is 
distinct from the structure of adult tissues. In addition, previous 
characterization of the early EB microenvironment indicated min- 
imal ECM deposition, largely devoid of collagen, during the first 4-7 
days of differentiation 44 ; thus, the principal structural attributes of 
stem cell aggregates are initially constituted by the cells and their 
cytoskeletal filaments, as well as by the intercellular adhesions (i.e. E- 
cadherin) supporting EB formation and maintaining multicellular 
organization. Consistent with our observations, a recent report 
describing a novel method to measure local forces within cellular 
aggregates attributed local stresses to myosin activity 45 . Moreover, 
the stiffness of embryonic tissue explants has been measured to be in 
the range of 10-100 Pa 14 , consistent with our reported data, thereby 
highlighting the unique physical traits of embryonic microenviron- 
ments. This study establishes a biophysical model of multicellular 
embryo-like tissues in which cytoskeletal tension is transmitted 
between cells via intercellular adhesions and is dynamically modu- 
lated during ESC differentiation and morphogenesis. 

In addition to the unique composition of embryonic-derived tis- 
sue microenvironments, pluripotent stem cells also exhibit distinct- 
ive biophysical characteristics, in contrast to their differentiated 
progeny. For example, hESCs exhibited a >3 fold increase in modu- 
lus upon differentiation toward chondrogenic lineages 46 and embry- 
onic tissues stiffen during gastrulation, leading to a 10-fold increase 
in mesoderm stiffness compared to endoderm 14 . Similarly, the data 
presented in our study collectively suggest that stiffness and viscosity 
increase concomitant with three-dimensional mesoderm differenti- 
ation and mesenchymal morphogenesis. The mechanism for 
increased stiffness upon differentiation has been suggested to result 
in part from cytoskeletal-mediated changes in structure 46 ' 47 . In addi- 
tion, recent studies have demonstrated the increasing stiffness of 
nuclear structures during mesenchymal stem cell differentiation 48 , 
with the biophysical characteristics attributed to lamin- A expression 
and chromatin structure 49 , which could account for the increasing 
contribution of factors independent of ROCK-mediated contractility 
at later stages of differentiation (Fig. 4d). Morphogenic events 
responsible for tissue-scale patterning, including local cell prolifera- 
tion 50 and EMT 33 , occur in regions of increased tension, indicating a 
role for biophysical regulation of stem cell differentiation and mor- 
phogenesis. In this study, perturbation of intracellular cytoskeletal 
elements during mesenchymal morphogenesis elucidated the rela- 
tive contribution to biophysical tissue properties; in ESC microen- 
vironments, the actomyosin cytoskeleton significantly contributed 
(>50% of modulus) to the microtissue stiffness during all stages of 
differentiation (Fig. 4a-c). Interestingly, early changes in stiffness 
and cytoskeletal organization have been used to successfully predict 
the differentiated phenotypes of adipose stem cells (ASCs) and 
MSCs, respectively 51,52 . Similarly, differences in cytoskeletal tension 
between basal and BMP4-treated EBs were detected as early as day 2 
of differentiation, suggesting that dynamic cytoskeletal changes 
occur prior to phenotypic differences commonly assessed by stand- 
ard molecular, biochemical or histological approaches. In addition, 
the assessment of biophysical microtissue characteristics provides a 
direct measurement of the dynamics of active processes mediating 
morphogenesis, thereby providing metrics to monitor processes such 
as cell sorting and patterning, analogous to similar events during 
embryogenesis 11,12 . Collectively, the data presented herein suggest 
that stem cell differentiation is reflected by concomitant changes in 
biophysical properties, thus establishing a novel, non-destructive 
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quantitative metric to assess microtissue morphogenesis and the 
homogeneity/heterogeneity of multicellular aggregate populations. 

Biomechanical regulation of stem cell fate has become a prom- 
inent research focus as an approach to engineer stem cell niches 23 . 
Manipulation of extracellular mechanical properties has been 
explored to direct stem cell differentiation, with changes in morpho- 
logy induced by modulation of cell spreading 37 or substrate stiffness 24 
linked to altered MSC differentiation trajectories. However, moni- 
toring and controlling microtissue biomechanics is particularly chal- 
lenging due to the limitations associated with manipulating or 
presenting three-dimensional biophysical cues 53 and an incomplete 
understanding of the physical mechanisms governing 3D stem cell 
differentiation and patterning. Large-scale bioinformatics tech- 
niques based on gene expression signatures are commonly used to 
define stem cell states 54 , but such approaches rely on data generated 
from a single, often destructive, analytical technique. The bio- 
mechanical assay demonstrated herein establishes an orthogonal, 
non-destructive metric for describing three-dimensional stem cell 
aggregate characteristics. Multivariate modeling (Fig. 5) established 
a quantitative relationship between gene expression and biomecha- 
nical properties of stem cell aggregates differentiated in various bio- 
chemical contexts, including those directed toward ectoderm (basal) 
and mesoderm lineages (BMP4), thereby providing proof-of-prin- 
ciple for the utility of biomechanical metrics to monitor stem cell fate 
by using non-destructive techniques to measure intact aggregate 
properties. Therefore, this study presents a novel combination of 
biophysical metrics to assess the 3D differentiation of ESCs, which 
highlights opportunities to study the biophysical characteristics of 
divergent differentiation trajectories, cell phenotypes, and morpho- 
genesis arising in stem cell-derived tissues. Together, the data pre- 
sented suggest that the phenotype of multicellular mesodermal 
microtissues can be described by both molecular (gene/protein) 
and/or biophysical signature 55 , which together define the morpho- 
genic state of three-dimensional pluripotent stem cell aggregates 
undergoing differentiation. Ultimately, the dynamic changes in 
embryonic tissue microenvironments illustrates important princi- 
ples for the design of strategies to direct differentiation and pattern- 
ing of complex stem cell-derived tissues amenable to modeling 
embryogenesis, screening pharmacological compounds, and devel- 
oping molecular therapeutics. 

Methods 

Embryonic stem cell expansion. Pluripotent murine ESCs (D3 cell line) were 
expanded on 0.1% gelatin coated polystyrene tissue culture treated dishes in 
Dulbecco's Modified Eagle's Medium (DMEM) containing 15% fetal bovine serum 
(FBS) and supplemented with penicillin (100 U/mL), streptomycin (100 mg), 
amphotericin (0.25 mg/mL), L-glutamine (2 raM), MEM non-essential amino acid 
solution (IX), 2-mercaptoethanol (0.1 mM), and leukemia inhibitory factor (LIF; 
10 3 U/mL). Media was exchanged every other day and ESCs were passaged prior to 
70% confluence. 

Embryoid body formation and culture. A single cell suspension of ESCs was 
obtained by treatment of monolayer cultures with 0.05% trypsin/0.53 mM EDTA 
solution. EBs were formed by centrifugation (200 rcf) of ESCs into 400 um diameter 
polydimethylsiloxane (PDMS) microwells {Aggrewell; 1M Stem Cell Technologies), 
with approximately 1000 cells per well. After 20 hours of micro well formation in 
serum containing ESC media without LIF, EBs were transferred to suspension culture 
(—2000 EBs per dish) and maintained on a rotary orbital shaker platform at 65 rpm, 
as described previously 26 . EBs were differentiated in basal serum-free media (N2B27) 
composed of DMEM/F12 supplemented with N2 {5 ug/ml insulin, 100 ug/ml 
transferrin, 6 ng/ml progesterone, 16 ug/ml putrescine, 30 nM sodium selenite) and 
50 ug/mL bovine serum albumin (BSA), combined 1 : 1 with Neurobasal IM media 
supplemented with B27 56 . Antibiotics and L-glutamine were supplemented at the 
same concentrations as the undifferentiated ESC media. 90% of EB media, including 
growth factor supplements, was exchanged every other day following gravity- induced 
sedimentation of EBs. Mesoderm induction was accomplished, as detailed in the 
results, by supplementation with 10 ng/mL BMP-4 (R&D Technologies) 57 . 
Perturbation of the cytoskeleton was accomplished by a 1-hour treatment with an 
inhibitor of the Rho- associated protein kinase (ROCK) pathway (Y27632) 32 , using the 
standard concentration established for ESC culture {10 uM) 58 or by supplementation 
of the media between days 4-7 of differentiation with Y27632, Jasplakinolide (Jas, 
50 nM) 34 , or Latrunculin B {LatB, 100 nM) 35 . In all cases, the concentration range for 



inhibitors was established based upon published measurements of ED 50 
concentrations 35 , as well as concentrations employed for extended culture durations 
(24+ hours) 34 . Moreover, concentrations were screened to ensure that small molecule 
treatments did not lead to significant loss of ESC viability and supported EB 
maintenance. DMSO (<0.1%) was included as a vehicle control. 

Parallel plate mechanical compression. The bulk mechanical properties of 
microtissues were measured using a micron-scale mechanical testing system 
(Microsquisher; CellScale), which calculates force via measurement of beam 
deflection in response to user-defined displacements (Supplementary Fig. Sla). All 
samples were tested in a PBS fluid bath (pH 7.4, containing 0.90 mM Ca 2+ and 
0.49 mM Mg 2+ ). The cantilever beams were composed of Tungsten (modulus — 
411 GPa), and beams of diameters from 76.2 um to 152.4 um were employed, 
depending on the stiffness and sensitivity required to measure different aggregates. 

To determine the bulk physical properties of EBs, a viscoelastic creep method was 
employed, which exerts a constant force while measuring the changes in displacement 
over time (Supplementary Fig. Sib). The magnitude of force was chosen as the 
average force corresponding to approximately 40% strain, determined based upon 
constant strain rate analysis of stress versus strain in n — 3 samples. The time to 
relaxation was determined to be less than two minutes, which was established based 
upon standard durations from published literature 46 ' 59 and from empirical testing, 
which established the steady-state phase based upon the constraint of less than 1% 
change in deformation over 1 second intervals. The creep parameter was calculated as 
the percent increase in displacement during constant force interval (Supplementary 
Fig. Sic). Relevant physical properties, including moduli (instantaneous modulus, E D ; 
relaxed modulus E int ), time constants (creep time constant, i a ; stress relaxation time 
constant, t c ), and apparent viscosity (p.) were calculated based upon a linear viscoe- 
lastic model of creep displacement (u) behavior over time, described in equations (1- 
3), where D D is the initial EB diameter, the stress (<r) is calculated as the constant force 
normalized to relaxed cross sectional EB area, and H(t) is a unit step function 59 . 



«(t) = 



2aD 0 



1+1 — —1 \H(f) 



E — —E 



(i) 



(2) 



(3) 



Quantitative real time PCR. RNA was extracted from EBs using the RNeasy Mini kit 
(Qiagen Inc, Valencia, CA) and cDNA was synthesized using the RT 2 First Strand Kit 
(SABiosciences, Frederick, MD) with 900 ng of total RNA. To assess gene expression, 
a custom RT 2 Profiler™ PCR Array (SABiosciences) was used in conjunction with a 
MylQ iCycler (Bio-rad) and iQ SYBR Green Supermix (Bio-rad). The array consists 
of 26 genes relevant to pluripotency and differentiation toward the three germ 
lineages, as well as 3 housekeeping genes (Actb, GAPDH, Hsp90abl) and PCR 
efficiency controls (Supplementary Table Si). Alternatively, quantitative real time 
PCR was performed using primers for the hyaluronan synthase isoforms (HAS-1, 
HAS-2, and HAS-3), which were designed using Beacon Designer software 25 . PCR 
data were normalized to the geometric average of the three housekeeping genes and 
fold changes in expression were calculated from normalized Ct values, relative to 
expression of ESCs 60 . The Genesis software package was used to generate heatmap 
visualizations of the gene expression data and to calculate the two-dimensional 
hierarchical clustering, based on Euclidean distance and average linkage clustering. 

Histology and immunohistochemistry. EBs sampled from rotary orbital culture 
were rinsed in PBS and fixed with 10% formalin (4% paraformaldehyde) for 45 
minutes with rotation at room temperature. Fixed EBs were embedded within 
Histogel (Richard Allen Scientific), and processed through a series of ethanol and 
xylene rinses prior to paraffin embedding. Paraffin embedded EBs were sectioned 
with a thickness of 5 um using a rotary microtome (Microtom HM310) and mounted 
on charged slides (SuperFrost Plus, Fisher Scientific). Sections were deparaffinized 
prior to staining via ethanol and xylene rinses and subsequently stained with 
hematoxylin and eosin, alcian blue (pH 2.5; nuclear fast red counterstain) or safranin- 
O (fast green counterstain). 

Alternatively, for immunohistochemistry, slides were blocked in 1.5% normal goat 
serum (NGS) for 1 hour at room temperature, incubated in mouse monoclonal anti- 
human smooth muscle actin (a-SMA) primary antibody (Dako; clone 1A4) for 1 
hour, and subsequently incubated in biotinylated donkey anti-mouse IgG (1 :400, 
Vector Labs). Slides were then incubated with avidin and biotin-based horseradish 
peroxidase, based on manufacturer instructions (Elite Vectastain Kit, Vector Labs) 
and developed using a 3, 3'-diaminobenzidine (DAB) substrate (Vector Labs). Nuclei 
were counterstained with hematoxylin. After staining, slides were cover slipped and 
imaged using a Nikon 80i upright microscope. Sections of mouse intestine, cartilage, 
and heart, were included as positive controls for alcian blue, safranin-O, and ot-SMA 
staining, respectively. 
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Multivariate partial least squares (PLS) regression modeling. A data matrix (MxN) 
was constructed with M culture conditions (Jas, LatB, etc), and N experimental 
parameters (Oct-4, Nanog, modulus, etc), where each column of the independent X 
(input) matrix represents a unique gene expression signal and the dependent 
variables in the Y (output) matrix represent mechanical properties. All data were 
mean-centered and scaled to unit variance and SIMCA-P (Umetrics) was employed 
to construct the PLSR model using a nonlinear iterative partial least squares 
(NIPALS) algorithm 6 '. In order to avoid overfitting the model, the number of 
significant principal components was determined based on both cross-validation 
(Q 2 Y) and permutation testing of each Y variable, with the background Q 2 Y less than 

0. 05 considered significant. Predictions were determined based on cross-validation, 
with standard error of predicted values determined via j ack-knifing, and the goodness 
of fit for the predicted data was calculated from coefficient of determination (R 2 ) 62 . 
Statistical significance testing was conducted by calculating the p-values from 
Pearson's r correlation coefficients of the predicted and observed data. Predictability 
was calculated from the root mean squared error of estimation (RMSEE) 63 . 

Statistical analysis. All experiments conducted with replicate data are represented as 
the mean of independent replicates + / — standard error, unless otherwise noted, with 
the experimental size (n) individually detailed in the corresponding figure legends. 
Prior to statistical analysis, all data were pre-processed using a box-cox power 
transform, in order to normalize data according to a Gaussian distribution. Statistical 
tests were conducted between groups and time points using one-way or two-way 
ANOVA, combined with either a post-hoc Tukey or Mann- Whitney U test for 
comparison of individual samples, depending on the results from Levene's equality of 
variances test, with p-values less than 0.05 considered significant. 
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